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Recent neutrino oscillation results have evoked renewed interest in sterile neutrino oscillation models. This paper 
reviews the data from MiniBooNE and short- baseline reactor antincutrino experiments within the context of 
sterile neutrinos. The results are incorporated into combined fits to test the viability of sterile neutrino oscillation 
models, which are later expanded to address matter effects. Finally, future experiments that can resolve the 
questions that have been raised are discussed. 



1. Sterile Neutrino Oscillation Formalism 

Sterile neutrinos refer to additional neutrino "flavor" (and mass) eigenstates which do not couple to the W ± 
and Z bosons. The additional mass eigenstates are assumed to include small admixtures of the three weakly 
interacting states (y e , v^, and v T ), so that they can be produced through mixing with the standard model 
neutrinos, and can therefore affect neutrino oscillations. 

In the case where only one such extra state is introduced to the standard, three-neutrino picture, the oscillation 
effects manifest as either v a disappearance or v a — > vp appearance, where a,/3 = e, n, or r. When the extra 

^ i mass eigenstate has mass 3> m^,mi, mi, the corresponding oscillation probabilities can be approximated 

by the two-neutrino oscillation probabilities, 

^ P^vl^vl) = 1 - sin 2 29 aa sin 2 (1.27 Am^L/E), (1) 

Oh 

■ and 
Qh 

Pfal^vfr*) = sin 2 26 afj sin 2 (1.27 Am^LfE), (2) 

for the case of disappearance and appearance, respectively. In the above expressions, the probability amplitudes 
sin 2 29 aa — 4|?7 a 4| 2 (l — |C/ Q 4| 2 ) and sin 2 29 a p — 4|J7 Q 4| 2 |(7 ( 94| 2 are functions which depend on the a and /3 flavor 
content of the extra mass eigenstate, |J7 Q 4| 2 and |J7/3 4 | 2 - 

In reality, more than one such state can exist. Indeed, models with two sterile neutrinos with to 5 ~ m 4 have 
C — ■ been extensively studied, specifically because of the possibility of observable CP violation effects which can arise 
through the 4-5 interference term which appears in the appearance 1 oscillation probability in this scenario pQ. 
The probability can be approximated for three-neutrino mixing as 

y—i 

^ P(t2^vp}a) = 4|C/ tt4 | 2 |C/ /}4 | 2 sin 2 (1.27ATO 2 lJ L/£;) 

+4|C/ a5 | 2 |^ 5 | 2 sin 2 (1.27Am 2 1 L/^) 

+8|C/ a5 ||(7 / 3 5 ||C/ a 4||^4|sin(1.27ATO^ 1 i/^)sin(1.27Am^ 1 L/S)cos(1.27ATO 2 4 i/^T^45)(3) 
where ^45 is a Dirac CP-violating phase, given by 

045 - arg(U* 5 U 05 U a4 U* 04 ). (4) 

2. Experimental Hints 

At present, the most significant evidence for sterile neutrino oscillations at Am 41 ~ 1 eV 2 and sin 2 29 ^ e ~ 
0.005 comes from the LSND experiment [2H1]. LSND was a tt + — > /i + decay-at-rest experiment which employed a 
liquid scintillator detector to identify D e events from possible P M oscillations by detecting the delayed coincidence 
of scintillation light emitted by a positron in the v e p —¥ e + n reaction and light emitted by the subsequent n 



Note that the disappearance probability in this scenario is insensitive to CP violation. 
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capture. In a search for P M — > v e oscillations at L/E ~ 0.4-1.4 m/MeV, LSND observed an excess of v e events 
above background prediction which was consistent with oscillations described by Eq. [2] at 3.8a. 

MiniBooNE, a pion decay-in-flight experiment, was designed to test the LSND evidence for oscillations using a 
higher-energy beam, and look for v e appearance at the same L/E as LSND. In contrast to LSND, MiniBooNE 
employed a higher energy neutrino beam and longer baseline, and relied on the Cherenkov signature of the e~ 
produced in a v e n — » e~p interaction to identify any v e resulting from oscillations of in the beam. The different 
neutrino energy, beam, and detector systematics, as well as the different event signatures and background at 
MiniBooNE, allowed for an independent test of the LSND result. Under the assumption of CPT conservation, 
which implies that Eq. [2] is identical for neutrinos and antineutrinos, the MiniBooNE search yielded no excess 
which could be consistent with the LSND signal [5]. Such excess would appear most noticeably in the 475- 
1250 MeV energy range in reconstructed electron neutrino energy. Instead, the observed v e data in that energy 
range were found consistent with background prediction. However, a 3. Oct excess of v e was observed at lower 
energy, below 475 MeV, corresponding to a significantly different L/E than that of the LSND excess. The 
lack of excess above 475 MeV in the MiniBooNE — > v e search resulted in an exclusion of the LSND 90% 
confidence level (CL) allowed parameters 2 at the 98% CL. At the same time, the low energy excess observed by 
MiniBooNE was found inconsistent with the single sterile neutrino oscillation scenario [BJ. 

Past attempts to simultaneously interpret the MiniBooNE neutrino and LSND antincutrino excesses included 
fits to CP-violating sterile neutrino oscillations 015], which are also discussed in Sec.[l] Other interpretations 
include heavy sterile neutrino decay [9"rUl| and other non-standard physics scenarios such as altered neutrino 
dispersion relations [12], extra dimensions [13], CPT- violation [14] . and non-standard neutrino interactions 
[15] 116] . At the same time, the possibility of mis-estimated background has been considered and ruled out by 
independent measurements at MiniBooNE [BJ, in particular for the case of NC single-photon 3 events, which 
make up the most dominant background to the MiniBooNE v e appearance search at low energy. However, the 
possibility of a new background, such as that contributed by anomaly- mediated single-photon production [17] . 
remains a viable interpretation for the MiniBooNE low energy excess. 

More recently, MiniBooNE performed a less-sensitive search for — > v e appearance [18] by inverting the 
neutrino beam polarity, resulting in a lower-intensity beam with similar energy as the beam. This search 
provided a direct test of the LSND signal interpreted as v e appearance with an oscillation probability given 
by Eq. [2] independent of any CPT or CP assumptions. Primarily due to low statistics, but also because of a 
small, ~ Ict excess observed in the 475-1250 MeV energy range, with an oscillation-like L/E dependence, the 
MiniBooNE — > v e search was found to favor oscillations over the null hypothesis at the 91% CL. At the same 
time, a low energy excess was also observed, though not as significant as that of the — > v e search. 

The above excesses are all observed in appearance channels, and if sterile neutrino oscillations were to be 
their underlying source, they would suggest observable v e and/or disappearance. Muon neutrino disap- 
pearance at Am 2 ~ 0.1 — 10 eV 2 is constrained to be less than 10% 4 by short-baseline experiments including 
MiniBooNE/SciBooNE [T5], CDHS [2U], and CCFR [21], which look for shape and normalization distortions in 
their respective event distributions as a function of neutrino energy, as well as MINOS [22], K2K and atmo- 
spheric neutrino experiments |23| . which look for normalization effects in their event distributions. Unless large 
v e disappearance is invoked, the disappearance constraints provided by the above experiments rule out the 
sterile neutrino oscillation interpretation of the LSND and MiniBooNE excesses. 

Electron neutrino disappearance, on the other hand, has been thought to be constrained to less than 10% by 
short-baseline reactor antincutrino experiments including Bugey [23] and CHOOZ [2S]. Recently, however, a re- 
analysis of past reactor antineutrino experiments, including measurements performed at ILL, Bugey, ROVNO, 
CHOOZ, and others (26], finds an overall normalization deficit observed by all those experiments, which is 
consistent with light sterile neutrino oscillations described by Eq. [IJ with Am 2 ~ 1 eV 2 and sin 2 29 ~ 0.1. The 
oscillation fits performed in [26) account for a recent re-evaluation of predicted reactor antineutrino flux spectra 
[27] . by which the new flux predictions are found to be ~3% higher than previously calculated, so that the 
past 9 e flux measurements performed by the above experiments are all now systematically 3% lower than the 
new predictions. The observed v e deficit corresponds to 2.14ct, and is referred to as the "reactor antineutrino 
anomaly." 

This newly-found possible hint for sterile neutrino oscillations, this time in the disappearance sector, 5 has 



2 Using Eq.pl and assuming CPT conservation. 

3 Cherenkov aetectors such as MiniBooNE cannot distinguish a single photon from a single electron. 

4 Approximate 90% CL upper limit to sin 2 20 MM . 

5 Note that a long-standing piece of evidence for much larger v a disappearance at very short baselines from the Gallex and SAGE 
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Figure 1: 90% and 99% CL alowed oscillation parameter regions obtained from fits to short-baseline antineutrino-only 
data. Left: Fit to v e disappearance data from CHOOZ and Bugey, accounting for the recent re-evaluation of predicted 
reactor i/ e fluxes [26]. Middle: Fit to i/ M — » P e appearance data from LSND and MiniBooNE. Also shown are the limits 
from the KARMEN — > i* e oscillation search. Right: Global fit to v e disappearance and — > S> e appearance data 
sets. The global antineutrino-only fit is consistent with short-baseline (3+1) oscillations described by Am^ ~ 1 eV 2 
and sin 2 26> M ~ 0.005 at > 3a CL. 



renewed the interest in global sterile neutrino oscillation fits [2"9"1 - I32"] . In Sec.|3j we present global oscillation fit 
results performed independently, 6 which confirm the findings in |32j . 



3. Review of Global Sterile Neutrino Oscillation Fits 

3.1. (3+1) Oscillations 

Under the single sterile neutrino oscillation scenario, also referred to as the (3+1) scenario, the MiniBooNE 
antineutrino and LSND excesses are indeed highly compatible with each other, as well as with all other short- 
baseline antineutrino data sets, including the KARMEN v„ — > v e appearance [33J and the Bugey and CHOOZ 
v e disappearance data sets. A combined analysis of the latter is also consistent with oscillations at the > 99% 
level. A "Parameter Goodness-of-fit" (PG) test [31] results in a 22% compatibility among all above antineutrino 
experiments. The level of compatibility among antineutrino data sets is illustrated qualitatively in Fig. [I] 

Of course, except for the purpose of making an empirical observation, namely the good compatibility of 
all antineutrino data sets under this simple oscillation assumption, one cannot ignore the constraints from 
the Vn disappearance experiments, as well as additional v e appearance constraints from NOMAD |35) and the 
MiniBooNE — > v e search. Under a CPT conserving (3+1) scenario, neutrino and antineutrino oscillations are 
identical, and therefore constraints from neutrino data sets are directly applicable to the global fits considered 
under this scenario. The level of constraint provided by the neutrino data sets is illustrated in Fig. [2j As 
suggested by the figure, it is impossible to reconcile all short-baseline results under a CPT-conserving (3+1) 
scenario. In this case, the compatibility of all short-baseline datasets returned by a PG test is 0.11%, essentially 
ruling out this oscillation hypothesis. 

3.2. (3-1-2) Oscillations 

As mentioned earlier, a (3+2) oscillation scenario is attractive because it offers the possibility of using CP 
violation to account for the observed differences in — what could be evidence of — short-baseline neutrino and 
antineutrino oscillations. Indeed, a combined analysis of all appearance data sets under the (3+2) scenario 
yields high compatibility [8] if one allows for large CP violation; furthermore, it predicts an excess of low energy 



calibration source experiments has been recently excluded by a comparison of KARMEN and LSND v e cross section measurements 
to theoretical predictions |28l . Those measurements, however, do not have sufficient sensitivity to address the reactor antineutrino 
anomaly. 

"From Ref. [B], with updates to include the updated reactor and MiniBooNE results, as noted. 



4 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 




10" 3 10" 2 10" 1 



sin 2 (2B^ e ) 

Figure 2: 90% and 99% CL alowed oscillation parameters (filled regions) and limits (lines) obtained from fits to short- 
baseline neutrino and antineutrino data. The solid lines correspond to 90% and 99% CL exclusion limits obtained from 
combined fits to null atmospheric disappearance and null short-baseline experiments. The filled region corresponds 
to the jointly-allowed LSND and MiniBoo NE antineutrino fit oscillation parameters, while the filled hashed region 
corresponds to the parameter region (marginally) allowed by MiniBooNE neutrino data. The figure is based on the 
analysis performed in [8], using an older, lower-statistics MiniBooNE antineutrino data set, and reactor v e disappearance 
data sets prior to the recent reactor flux re-analysis. Note that the updates do not significantly alter the conclusions 
drawn from this picture. 



events at MiniBooNE in both neutrino and antineutrino running, which can account for the observed excesses, 
as illustrated in Fig. [3} However, when disappearance data sets all included in the fits, the compatibility reduces 
dramatically, and one finds large tension among neutrino and antineutrino data sets. 



4. Expanding the Formalism to Account for Matter-like Effects 

The present (3+1) and (3+2) sterile neutrino oscillation fits reveal that neutrino and antineutrino oscilla- 
tion signatures which may be responsible for the LSND, MiniBooNE, and reactor antineutrino anomalies are 
incompatible with other, primarily neutrino, null results. The neutrino and antineutrino data sets show differ- 
ences which cannot be explained away by just CP violation. That observation, as well as a (now essentially 
resolved) hint from the MINOS disappearance search suggesting CPT-violating differences in neutrino and 
antineutrino atmospheric oscillations [361 137) , has prompted the theoretical community to consider even less 
standard scenarios, beyond the "reference picture" of CPT-conserving 3+N oscillations [58] . An example of 
such scenarios, involving sterile neutrinos and non-standard matter-like effects |39j . is discussed below. 

4.1. Extended Oscillation Formalism 

The hypothesis being investigated in this section is one where the appearance signals from MiniBooNE 

neutrino and antineutrino running and LSND are due to a (3+1) oscillation scenario, where W experience a 
matter-like potential V s = ±A S . The resulting appearance oscillation probability under this scenario is given 
by 

p(v^ v}) = sin 2 20 M sin 2 (l.27Am 2 M L/E), (5) 

where sin 2 26m and Am 2 M are the effective mixing parameters due to the presence of V s + 1 0. The oscillation 
probability has a form identical to that of a (3+1) scenario, except that the effective oscillation parameters 
are dependent both on the underlying "vacuum" (3+1) oscillation parameters and the neutrino energy and V s 
potential, according to 



Amfvf = Aml 1 +2EV S , 



(6) 
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Figure 3: The MiniBooNE v e (left) and £> e (right) event distributions versus reconstructed neutrino energy, from [8]. The 
null (no-oscillations) predictions are shown in green, with gray systematic error bands. The observed data are shown 
in black points, with statistical error bars. The signal predictions from a joint (3+2) CP-violating fit to all appearance 
(dashed red) and all appearance and disappearance (solid red) data sets are overlaid. The (3+2) CP-violating fit to all 
data sets fails to explain the observed excess signatures at MiniBooNE. 



and 



sin 2 9 M 



((Amj, - 2EV S ) 2 + 8EV s Am 2 41 \U s4 \ 2 ) 



(2EV S - Am^l - 2|[/ s4 | 2 ) + ^ (Am\ x - 2EV S ) 2 + 8£KA^|^f) 2 



(7) 



The dependence on EV S , where V s — ±A S , suggests different and also energy-dependent effective oscillation 
parameters for neutrinos versus antineutrinos. More specifically, for a set of vacuum oscillation parameters 
Am 2 !, It^l 2 , and |£/ e 4| 2 , a non-zero A s will lead to opposing deviations from the vacuum oscillation parameters 
for neutrinos versus antineutrinos, as well as resonance effects which shift to lower neutrino and antineutrino 
energies as A s increases. 

The following subsection presents results from a joint fit to this scenario performed using the MiniBooNE 
and LSND appearance data sets, assuming — > v e and — >■ v e oscillations described by Eq. [5] During the 
fit, the vacuum oscillation parameters, Am 2 1; |[/e,^4| 2 , and A s , are varied within 0.01-100 eV 7 , 0-0.05, and 
10~ 13 -10~ 9 eV, respectively. 7 



4.2. Fit Results 



Relative to a (3+1) fit, the compatibility among MiniBooNE and LSND data sets when one allows for non- 
zero A s increases from 2.3% to 17.4%. The allowed oscillation parameters are shown in Fig. [31 The returned 
best-fit vaccum oscillation parameters corespond to Am^ = 0.47 eV 2 and sin 2 2# Me = 4|[/ e 4p|[/ M 4| 2 = 0.01, 
and the best-fit A s value corresponds to 2xl0 -10 eV. Note that the best-fit A s value is significantly larger than 
A e = \/2Gi?rj e ^ 10~ 13 eV for standard matter effects. 

The effective mixing parameters corresponding to the best-fit vacuum oscillation parameters and A s are shown 
as a function of neutrino (and antineutrino) energy in Fig. [5] The effective parameters at lower energies, below 10 
MeV, reduce to the vacuum (3+1) oscillation parameters, and could likely accommodate the reactor antineutrino 
anomaly (3+1) oscillation interpretation. At higher energy, however, the effective mixing parameters suggest 
observable effects in the MINOS, K2K, IceCube, and atmospheric antineutrino measurements. Those effects 
are currently being investigated [39] . 



The choice of A a range is driven by short-baseline data. 
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Figure 4: The 90% and 99% CL allowed vacuum oscillation parameters Am 2 ! and sin 2 29^ e = 4|[/ e 4| 2 |t^j4| 2 an d A s 
potential obtained from a joint (3+1) with matter-like effects fit to the MiniBooNE and LSND data sets. The star 
indicates the best-fit parameters. 
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Figure 5: Energy dependence of the effective mixing parameters corresponding to the best-fit parameters in Fig. [5] The 
gray flat lines correspond to the best-fit vacuum oscillation parameters, while the solid (dashed) black lines correspond 
to the effective best-fit oscillation parameters for neutrinos (antineutrinos) . The vertical lines approximate the mean 
neutrino energy in the MiniBooNE and LSND antineutrino (dashed) and neutrino (solid) appearance searches. At low 
energies, the effective oscillation parameters reduce to the vacuum ones, and no differences in neutrino versus antineutrino 
oscillations are expected; at higher energies, one sees increasing deviations from the vacuum case as well as increasingly 
larger differences between neutrinos and antineutrinos. The plot suggests observable effects at high-energy (> lGeV) 
antineutrino data sets. 



5. Future Tests 

Experiments which aim to address the short-baseline anomalies discussed in Sec. [2] include both appearance 
and disappearance experiments. Newly proposed disappearance experiments include reactor antineutrino dis- 
appearance searches at very short baselines of the order of 10 m or less, where the L/E dependence of possible 
oscillation effects can be directly observed rather than inferred due to an overall normalization deficit, as is 
presently the experimental case. Such experiments include NUCIFER 40 and SCRAAM [31]. Similar searches 
are proposed for very short baseline v e and z/ e disappearance using MegaCurie-strength radioactive sources 
placed either inside or in close proximity to a neutrino detector, such that oscillation "wiggles" can be fully 
observed within the detector spatial dimensions. An example of such experiment is discussed in [42) . 

Of course, appearance-style short-baseline oscillation experiments are also being proposed, where particular 
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effort is being made to focus on experiments with the highest possibie (ideally, 5a) sensitivity to — > v & and/or 
v^, — > v e appearance. Those proposals include short-baseline oscillation searches at OscSNS, BooNE, Gd-doped 
Super-K running with a decay-at-rest neutrino source, and coherent scattering experiments with a decay-at-rest 
neutrino source [HHIS]- 

At the same time, several presently running and planned experiments aim to shed more light onto the 
picture. Specifically, MiniBooNE antineutrino running has been ongoing, aiming to nearly double the present 
antineutrino statistics to a sample corresponding to 1.5xl0 20 POT by mid 2012. The higher-statistics search 
would still be statistics- rather than systematics- limited, but could nevertheless either strengthen the significance 
of the present antineutrino result to ~ 3er, or reduce it to the < la level. 

A more definitive test can be performed by MicroBooNE [47) . which is an experiment under construction 
which will run in the same neutrino beamline as MiniBooNE, at roughly the same baseline. In contrast to 
MiniBooNE, MicroBooNE uses a 60-ton (fiducial mass) liquid argon time projection chamber (LArTPC) de- 
tector, which has a significantly higher detection and PID efficiency than MiniBooNE. MicroBooNE presently 
plans to run in a neutrino beam, so that it can perform a definitive test of the MiniBooNE low energy excess. 
Specifically, MicroBooNE's single-electron versus single-photon high differentiation capability grants a > 5a 
sensitivity to a single-electron interpretation of the low energy excess events, and a > 4er sensitivity to a single- 
photon interpretation. MicroBooNE is also in position to perform a — » v e appearance search; however, its 
relatively small size limits the statistical sensitivity of such search to a level which would only be comparable 
to the MiniBooNE sensitivity. A proposed future MicroBooNE run with a second, larger LArTPC detector 
48J, however, could provide a > 5ct sensitivity in both neutrino running, as well as antineutrino running. A 
similar oscillation search involving two LArTPC detectors at the CERN-PS beam has also been proposed by 
an independent group [15] , 

6. Summary 

Over the last five years, a number of hints have surfaced which are in support of the LSND evidence for sterile 
neutrino oscillations. Those hints come from — > v e and — > v e appearance measurements at MiniBooNE, 
and a combined analysis of past D e disappearance measurements performed at reactors at short baselines. At the 
same time, however, constraints from disappearance measurements by MINOS, K2K, MiniBooNE, and other 
experiments which show no evidence for such oscillations, place stringent limits to those models. Over the past 
five years, global fit results have revealed increasing tension between neutrino and antineutrino data sets when 
tested under the sterile neutrino oscillation hypothesis. The apparent differences are such that they cannot 
be accommodated by any CP violation allowed within the model. Non-standard interactions and matter-like 
effects, such as of the type discussed in Sec. [5j generally induce effects which look like CPT violation, and 
therefore show promise in being able to simultaneously explain the LSND, MiniBooNE, and reactor anomalies. 
Those effects have received increased interest within the community over the past year, and should be further 
developed and tested. An example of such model has been presented in this paper. While the model yields a 
satisfactory fit to the data, it requires a significantly large "matter-like" potential V s = ±A S , experienced only 
by sterile neutrinos, which invites theoretical interpretation. At the same time, it suggests observable effects in 
higher-energy neutrino and antineutrino experiments. Therefore, the data provided by MINOS, IceCube and 
T2K will be crucial in addressing the viability of the model. On the other hand, high-precision, very short 
baseline experiments at reactors should be able to directly test the sterile neutrino oscillation interpretation for 
the reactor anomaly. Several proposals for new appearance-style experiments are also being developed, with 
the goal of a definitive, > 5a test of the LSND and MiniBooNE results. 
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